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Local structure characterization with the bond-orientational order parameters q±, q$,. . . introduced by 
Steinhardt et al. has become a standard tool in condensed matter physics, with applications including glass, 
jamming, melting or crystallization transitions and cluster formation. Here we discuss two fundamental flaws 
in the definition of these parameters that significantly affect their interpretation for studies of disordered 
systems, and offer a remedy. First, the definition of the bond-orientational order parameters considers the 
geometrical arrangement of a set of neighboring spheres NN(p) around a given central particle p; we show 
that procedure to select the spheres constituting neighborhood NN(p) can have greater influence on both 
the numerical values and qualitative trend of qi than a change of the physical parameters, such as packing 
fraction. Second, the discrete nature of neighborhood implies that NN(p) is not a continuous function 
of the particle coordinates; this discontinuity, inherited by qi, leads to a lack of robustness of the qi as 
structure metrics. Both issues can be avoided by a morphometric approach leading to the robust Minkowski 
structure metrics q[. These q[ are mathematically equivalent to the recently introduced Minkowski tensors 
[Europhys. Lett. 90, 34001 (2010); Phys. Rev. E. 85, 030301 (2012)]; they are of a similar mathematical 
form as the conventional bond-orientational order parameters. 
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In 1983 Steinhardt et alX proposed the family of lo- 
cal qi and global Qi bond-orientational order (BOO) pa- 
rameters as a three dimensional generalization of the i/jq 
hexatic order parameter in two dimensions^. Bond ori- 
entation analysis has become the most commonly used 
tool for the identification of different crystalline phases 
and clusters, notably fee, hep and bco 3 - - — or icosahe- 
dral nucleiiSr— . They are also used to study melting 
transition a 10 i 13 i 14 and interfaces in colloidal fluids and 
crystals^. For the study of glasses and super-cooled 
fluids g6 and Qe have become the most prominent or- 
der parameter when searching for glass transitions^— 
and crystalline cluste ra 4 - 8 ' 11 ! 20 " — . While qi is defined as 
a local parameter for each particle, other studies have 
used global averages of bond angles (Qi) to detect single- 
crystalline order across the entire sample 2 ^—. 

The BOO parameters qi and Qi are defined as struc- 
ture metrics for ensembles of N spherical particles. For 
a given sphere a one assigns a set of nearest neighbors 
(NN) spheres NN(a). The number of NN assigned to a 
is n(a) — |NN(a)|. Any two spheres a and b are said 
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to be connected by a bond if they are neighbors, i.e. if 
a G NN(6)2£. The set of all bonds is called the bond net- 
work. The idea of bond orientation analysis is to derive 
scalar metrics from the information of the bond network 
(i.e. the set of bond vectors). The precise definition of 
the bond network is therefore crucial. 

Other structure metrics are defined in a similar way, 
with the only difference that the geometric interpre- 
tation of the bond network is changed; e.g. centro- 
symmetry metrics 2 ^ or Edwards configurational tensors 2 ^ 
and fec/hep-order metrics 2 ^, or the number of bonds as 
the most simple topological characteristic* 3 -^ 

For a sphere a the set of unit vectors n a b point from a 
to the spheres b € NN(o) in the neighborhood of a. Each 
vector n a f, is characterized by its angles in spherical coor- 
dinates 9 a b and ip a b on the unit sphere. Following Stein- 
hardt et aZ»i, the local BOO qi(a) of weight I assigned to 
sphere a is defined as 
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where Yi m are spherical harmonics (see e.g. appendix 
in2i). This formula can be interpreted as the lowest- 



order rotation-invariant (that is, independent of the co- 
ordinate system in which Q ao and f a b are measured) of 
the Z-th-moment in a multipole expansion of the bond 
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vector distribution Pbond(n) on a unit sphere. Higher- 
order invariants, often termed wi, are defined in a similar 
wayi^— . 

The existence of spheres with values of and q§ 
close to those of an ideal ordered structure has been in- 
terpreted as evidence of ordered clusters (see Tab. HJ. 
The local structure metrics qi have been used to iden- 
tify fee, hep, bec or icosahedral structures in condensed 
matter and plasma physics (e. g. in colloidal particle 
systems*, random sphere packing s 23 ! 3 or plasmas^ 5 .) by 
analyzing histograms over the ((74, qe)-p\a,ne or combi- 
nations of similar order parameters^. Frequently, his- 
tograms of one order parameter only, namely q§ , are used 
to qualitatively compare disorder in particulate matter 
systems ^ 20 ' 36 ! 37 . Our previous worb 3 ^ has raised the cau- 
tion that local configurations can exist that are clearly 
non-crystalline but have the same values of q§ &s hep 
or fee environments. Several authors have defined bond 
order functions 3 ^ closely related to the qi for the identi- 
fication of crystalline cluster s 11 ! 15 ' 21 ! 40 . 

As a different application from the identification of 
locally crystalline domains, it has been proposed to 
use averages (qi) over all spheres to quantify the de- 
gree of order of a configuration. Averages (gg) have 
been analyzed (as function of some control parameter 
such as temperature, pressure, strain, or packing frac- 
tion) for random sphere packings 2 ^, granular packing 
experiments 4 ^, model fluids^ 2 -, molecular dynamics simu- 
lations of water— or polymer melts^ 4 -. This use of (qi) to 
quantify the overall degree of order implies a monotonous 
relationship between the value of qi and the degree of or- 
der. In contrast to the identification of individual crys- 
talline cells as those with qi the same as for the crystalline 
reference cell gf ryst , one now assumes that larger values 
of A := \qi~ Q; Clyst | correspond to "larger" deviations from 
the crystalline configuration, even for clearly acrystalline 
local configurations with large values of A. The validity 
of this assumption is difficult to assert, in the absence of 
an independent definition of the degree of the "deviation 
from crystalline structure". (Note also the obvious prob- 
lem for the case of monodisperse hard spheres, where two 
distinct crystal reference states, fee and hep, exist which 
however have different values of qi.) Nevertheless, <jg has 
been used to quantify order in disordered packings, un- 
der the assumption that higher values of q§ correspond 
to higher degree of order—. Unless the system represents 
a small perturbation of one specific crystalline state, this 
use of qe is, in our opinion, not justified. q$ is n ot a 
suitable order metric to compare the degree of order of 
disordered configurations that are far away from a crys- 
talline reference state. We use the term structure metric 
to emphasize that a priori qi does not quantify order in 
disordered systems. 

We here demonstrate a further aspect, distinct to those 
described above, that should be taken into account when 
interpreting qi data for disordered systems, namely a very 
significant dependence of the qi values on details of the 
definition of the bond network: changes of the NN defi- 
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TABLE I. Values of qi in perfectly symmetric configurations. 
For these highly symmetric cases (fee, hep, icosahedron,sc), 
the definitions of neighborhood discussed in this article all 
yield the same crystallographic neighbors, and hence values 
of qi (assuming infinite precision for the point coordinates 
such that the Delaunay diagram has edges to all nearest crys- 
tallographic neighbors). Spheres in bee configuration have 8 
nearest neighbors at distance a, where a is the particle di- 
ameter, and 6 second nearest neighbors at distance y/2a and 
have 14 Delaunay neighbors. 

nition do not only affect the absolute values (which are 
of great importance, as the comparison to the crystalline 
reference values is in terms of these absolute values) but 
they can also affect functional trends. This observation 
highlights the problem in the interpretation of anomalies 
of the BOO parameters (that is, local extrema as function 
of some thermodynamic parameter) as being connected 
to thermodynamic anomalie n 42 ' 43 ; see also the discussion 
of the anomalies of water— in terms of a parameter sim- 
ilar to the BOO parameters. Rather than being a mere 
inconvenience, the dependence on the details of the bond 
network definition is of direct relevance to the physical 
interpretation. 



AMBIGUITY OF THE NEIGHBORHOOD DEFINITION 
AND ITS EFFECT ON qi 

The choice of a set of nearest neighbors - at the heart 
of bond orientation analysis - is not unique (see Fig. [IJ. 
Steinhardt et al. proposed to use "some suitable set" of 
bonds for the computation of qi\ they used a definition 
based on a cutoff radius of 1.2a, where a is the parti- 
cle diameter—. That is, each sphere that is closer to a 
given sphere a than a cutoff radius r c is assigned as a 
NN of sphere a. Neighborhood definitions based on cut- 
off radii are widely used, e.g. with cutoff radii 1.2a and 
\Ao±±^£2£L£L or with the value of the cutoff radius 
determined by the the first minimum of the two-point 
correlation function g(r) & 14 i 15 i 25 i 48 . 

Alternatively, the Delaunay graph of the particle 
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FIG. 1. (color online) Widely used NN definitions: a) Voronoi 
diagram (red) and its dual, the Delaunay graph (blue) b) 
Delaunay definition of nearest neighbors (NN): the Delaunay 
neighbors of the red sphere are highlighted in green, c) NN 
definition with cutoff radius r c d) rif closest NN, here rif — 6 



centers^ ^ is used to define NN^SSaSMiiSi. In this 
parameter-free method, every sphere which is connected 
to a sphere a by a Delaunay edge is considered a NN 
of a. A rarely used definition is to assign a fixed num- 
ber rif of NN to each particle n{a) = rifi^S. In three 
dimensions, the rif = 12 other spheres closest to the cen- 
tral sphere are chosen as neighbors. The difference be- 
tween these definitions is illustrated in Fig. [T] Note that 
while the definitions via cutoff radius and via the Delau- 
nay graph, panels (b) and (c) in Fig. [TJ are symmetric, 
i.e. b € NN(o) ««6 NN(fc), the definition of neighbor- 
hood as the nearest rif spheres is not, see Fig.Q](d). The 
definitions of NN discussed so far will be called bond net- 
work neighborhoods in the following; in this picture, each 
nearest neighbor is equivalent to the other neighbors. By 
contrast, we use the term morphometric neighborhood if 
the neighborhood relation is additionally weighted with 
geometrical features. 

A principal weakness of structure metrics based on 
bond network neighborhoods is their lack of robustness: 
Small changes of particle positions can delete or add en- 
tries in the set of neighbors. This discontinuity w. r. t. the 
particle positions is inherited by the structure metrics de- 
fined via bond network neighborhoods. Small changes in 
the particle coordinates can lead to large changes in the 
structure metrics, which is undesirable. 

We demonstrate the very strong effect of the NN defini- 
tion on the BOO parameter qe by the example of a super- 
cooled fluid. Using non-equilibrium molecular dynamics 
(MD) simulations^^— , super-cooled configurations are 
generated that represent entirely disordered states with 
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FIG. 2. (color online) Average local bond order parameter 
(qe) in the super-cooled HS fluid with several definitions of the 
nearest neighbors: orange squares: r c = 1.2a, green bullets: 
r c = 1.4<t blue crosses: Delaunay definition and black stars: 
n f = 12. Turquois triangles based on (q' 6 ) defined in Eq. ([2}. 



densities larger than the fluid-crystal coexistence density 
of hard spheres (HS) of </> 0.494^. 

Figure [2] shows the average local BOO (qe) for four 
different choices of bond network neighborhood defini- 
tion. To distinguish between the different definitions 
of neighborhood discussed above, we use the symbols 
q r 6 " , q® and q§ { . First, the absolute values of <7g c_1 ' 2cr , 
%"~ 1A(J i % f ~ 12 and ^ip differ significantly, which is im- 
portant when comparing these values to that of a specific 
crystalline phase such as fee. Second, and of greater con- 
cern for the use of q§ as a structure metric, the behavior 
of q^ 1 - 2 *, q r 6 c=1Aa , ?6 f=12 and g° is qualitatively dif- 
ferent as a function of the packing fraction </>. For exam- 
ple (<?g c_1 ' 2c r )(</>) shows a slight negative trend without 
pronounced extrema, whilst (<?g c_1 ' 4<T }(</>) increases for 



4> < 0.56 and decreases above. {q§ t=l2 )((j>) and {q^){(f>) 
show a maximum at slightly different positions with a sig- 
nificantly different absolute value. Each of these trends 
is specific to the neighborhood definition. These discrep- 
ancies raise a caution flag about the use of gg as a local 
structure metric in disordered systems. This is in accor- 
dance to several reported difficulties in the application of 
gg in ordered and disordered system o 45 ' 56 ' 57 . The choice 
of the NN definition has a dominant effect on the values 
and on the functional trend of {qe)(4>) that conceals the 
behavior due to genuine structural changes induced by 
the physics of the system. Results for gg obtained by 
different studies are not only difficult to compare quan- 
titatively, but also the qualitative behavior may be mis- 
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FIG. 3. (color online) (a) Average number of nearest neigh- 
bors identified by the different definitions of neighborhood, 
for the same data as shown in Fig. [2 Difference (qg) - (q' 6 ) 
between qe values for different definitions of the bond net- 
work neighborhood X = {r c = 1.2a, r c = 1.4cr,D,nf = 12}. 
(b) Comparison of the functional trend of these data to 
— {n)((/>) demonstrates the strong negative correlation of the 
value of q§ with the number of NN spheres n identified by the 
specific neighborhood definition. 



leading. 

The behavior of (qe) can be rationalized by considering 
the average number of nearest neighbor spheres (n)(c/>) 
identified by the different neighborhood definitions. 

Figure [3] (a) shows (q 6 ) — (q' 6 ) as function of </>. q' 6 is a 
structure metric based on morphometric neighborhood, 
which is discussed in detail in the following section. Fig- 
ure [3] (b) shows — (n) . The very close correlation of the 
curves with the same bond network neighborhood defi- 
nition is evident. This result asserts that (q' e ) captures 
physical structure properties, while various variants of 
(cfe) are predominantly indicative of the typical number 
of NN spheres (n) identified by the respective NN defini- 
tions. 

Figure |4] further corroborates this observation by the 
analysis of (gg f- ™) as a function of n for the super-cooled 
hard sphere fluid at <f> = 0.600. The average (qQ t ~ n )(n) 
systematically decreases with higher prescribed numbers 
m of NN. This effect is further amplified for large m > 
12, when spheres in the second coordination shell are 
also identified as neighbors. The stronger decrease in 
qe when encountering the second coordination shell also 
explains why (q®) generally has lower values compared 
to the other neighborhood definitions, since } « 14. 



MINKOWSKI STRUCTURE METRIC BY VORONOI-CELL 
WEIGHTING 

This section introduces the Minkowski structure met- 
rics q[ that we already alluded to above. The Minkowski 
structure metrics (MSM) are obtained by an adaption 
of the conventional BOO parameters. The MSM differ 
from the conventional qi, Eq. (JTJ) , by the fact that the 
contribution of each neighbor to the structure metric is 
weighted by an associated relative area factor A(f)/A. 



FIG. 4. (color online) (a) Mean (q^^™) as a function of the 
fixed number n of neighbors assigned to each sphere. The 
squares are data of a super-cooled fluid with 4> = 0.6 and the 
red solid line of the ideal gas ({qg f ~ n ) oc n -1//2 , see Ref.— ). 
The dotted lines are fits for the first coordination shell for 
n < 12 and first and second shell n > 12. The first shell ex- 
ponent is —0.24 and the second shell exponent is —1.48. (b) Il- 
lustration for the discontinuity of the topology of the Voronoi 
diagram as function of center point coordinates: An infinites- 
imal particle displacement can destroy or create Voronoi cell 
facets (and hence bonds in the neighborhood definition based 
on the Delaunay graph). 



In this factor, A(f) is the surface area of the Voronoi 
cell facet / that separates two neighboring spheres, and 
A = Y^f£F( a ) A(f) is the total surface area of the Voronoi 
cell boundary J- (a) of sphere a. This simple change leads 
to robust, continuous and parameter-free structure met- 
rics q[ that avoid the shortcomings of the conventional qi 
discussed above. 
We define 
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where 8f and iff are the spherical angles of the outer 
normal vector n/ of facet /. Note that the direction of 
this vector coincides with the bond vector that is used in 
conventional bond orientation analysis (see Fig.[l|. 

Because of the weighting of each bond by its corre- 
sponding Voronoi facet area A(f)/A, these newly con- 
structed structure metrics q[ are continuous functions of 
the spheres' center point coordinates, and hence robust. 
Furthermore, this geometrical neighborhood is symmet- 
ric and parameter-free. 

The definition of q[ results naturally from a multipole 
expansion in spherical harmonics of the Voronoi cell sur- 
face normal distribution function 



p(n) = j.J2S(n(f)-n)A(f) 



(3) 



on the unit sphere: /°( n ) = p[@i ^P) = 
EtoEL-iWmf^), where q[ m evaluates to 
Y, f er(a)( A (f)/ A ) Y * m ( 9 f,Vf); the star denoting com- 
plex conjugation. 
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By contrast, the Zth-moment of the distribution p(n) 
in Cartesian coordinates is 
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where <X> denotes the tensor product. The moment ten- 
sors are special types of Minkowski tensor a 53 ' 59 . 
These versatile shape metrics are intensively studied in 
the field of integral geometry^ and successfully applied 
to analyze structure in jammed bead pack a 61 ' 62 , bi-phasic 
assemblie a 63 i 64 , foams and other cellular structure a 59 ' 65 . 
There is a one-to-one correspondence between this class 
of Minkowski tensors and the multipole expansion of 
the surface normal vector distribution p(n) of a convex 
Voronoi polytope T(a) 66 ^ 7 . 

For ideal crystals where all Voronoi facets have equal 
size, the values of the BOO qi and of the MSM q[ are the 
same; These symmetries are fee, hep, the icosahedron and 
sc (simple cubic). In the case of bec, where Voronoi cells 
have in total 14 facets, of which 8 correspond to closest 
neighbors and 6 to neighbors in the second shell, qi differ 
from q[ (see also Table |TJ) . 



GEOMETRIC INTERPRETATION OF THE MINKOWSKI 
STRUCTURE METRICS, IN PARTICULAR OF q 2 

For the use of both BOO parameters and MSM, an im- 
portant issue is the choice of the weights I that are con- 
sidered. Many studies restrict themselves to only q§, pos- 
sibly supplemented by q^ and the associated higher-order 
invariants and wq. This is likely to be motivated by gg 
being the apparent generalization of the two-dimensional 
hexatic order parameter if>e- The relation between the 
I = 6 structure metrics and ordering, however, is not as 
direct in 3D as it is in 2D: q§ is maximized by icosahedral 
bond order, which is incompatible with translational or- 
der. The perception that large values of certain structure 
metrics, in particular q§, are intrinsically connected with 
crystallization is therefore deceiving, and it is useful to 
discuss the relevance of the individual weights to physical 
problems. 

In all cases, q' is trivially one while q[ trivially van- 
ishes, due to the so-called envelope theorems of Mueller— 
(note, this does not apply to qi). Thus, the first weight 
that captures pertinent information about a disordered 
system is I = 2; for hep and fee crystals q' 2 vanishes. 
The invariants q^ and q' 5 (and odd weights in general) 
vanish in configurations symmetric under inversion, but 
capture deviations from this symmetry (see tab. [I]). Hence 
they might be robust candidates for defect detection like 
centro-symmetry metrics 2 ^ or to separate hep from fee, 
since the hep Voronoi cell is not inversion symmetric (see 
tab. H| , while fee is inversion symmetric {mim) with re- 
spect to the sphere centers. Including Steinhardt et al.'s 
original paper— we are not aware of any applications of 
odd weights I. The lowest weight to discriminate a sphere 
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FIG. 5. (color online) Minkowski structure met- 

rics q' 2 , q'i and q' 6 for equilibrium hard spheres (Monte 
Carlo, MC)^ simulations, jammed Lubachevsky-Stillinger 
(jLS)£i, non-equilibrium unjammed Lubachevsky-Stillinger 
(uLS)2. and non-equilibrium Matsumoto algorithm (MA- 
MD) simulations^ (see text). 0l' 2 is the anisotropy index, 
i.e. the ratio of the smallest and the largest eigenvalue of the 
Minkowski tensor W°' 2 ; see Eq. (QJ and Ref.— . 



from a cube is / = 4 and thus plays an important role 
in ordered materials. The cubic-symmetry fee, bec, and 
simple cubic lattices all have non- vanishing q' 4 values (for 
the conventional BOO parameters though, great care is 
needed for the bond definition, as different sets of NN for 
bec reveals a dramatic change on conventional 54). gg is 
the first non-vanishing weight for icosahedral symmetry 
(and maximum for the icosahedron). Note that the qe 
values for fee can be matched by deformed icosahedral 
bonds. 

While in ordered states, the qi are easily interpreted, in 
disordered states the lack of a well-defined reference state 
renders the interpretation more difficult. Fig. [5] shows 
(q' 2 ), (q'4) and (<7 6 ) of hard-sphere systems in a wide range 
of packing fractions. The plot includes data from Monte 
Carlo simulations of the thermal equilibrium fluid/solid^ 3 - 
(MC), from fully disordered and partially crystalline 
jammed Lubachevsky-Stillinger (jLS ) 61 ' 69 , and also from 
unjammed non-equilibrium simulations (uLS) from LS 
simulations before jamming^ and the data from Fig. [2] 
(MA-MD)^ 2 .. 

Empirically, we find that disordered cells virtually al- 
ways have finite q' 2 values; for order (cubic-symmetry or 
close packed), q' 2 vanishes. Therefore, distributions of q' 2 
in a partially ordered system are bimodal, which is conve- 
nient for the separation of both phases. Conversely, if the 
abundance of small values q 2 « in a sample vanishes, 
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one can conclude that it is fully disordered. The infor- 
mation contained in the lowest weight q' 2 is also captured 
in the anisotropy index J 8 1 ' derived from Minkowski 
tensors^. The observation that q 2 vanishes for ordered 
configurations corresponds to the fact that j3{ =1, and 
q' 2 > corresponds to /3°' 2 < 1 (cf. Refs.^&Si). 

Both structure metrics, q' 2 and /3°' 2 , capture well the 
different features in local structure of hard-sphere sys- 
tems (Fig. [5l panels (a) and (c)). The thermodynamic 
phase transition from the fluid to the solid (fee) phase 
at packing fractions around (f> w 0.49 is clearly visi- 
ble. Furthermore, jammed sphere packs are well dis- 
tinguished from the equilibrium configurations. Start- 
ing from the equilibrium and avoiding crystallization, 
the non-equilibrium MA-MD protocol continues the fluid 
branch into a super-cooled fluid regime. The uLS proto- 
col generates further non-equilibrium states with larger 
q' 2 , up to jammed configurations. In both diagrams (a) 
and (c), the non-equilibrium fluid states are found above 
the linear extrapolation of the equilibrium fluid branch, 
while the ordered phase is below. The diagram (b), show- 
ing (q' 6 ), reproduces (though "upside down") quite well 
the qualitative features obtained from {q' 2 } or {1 — j3{ ). 
The agreement of these two plots, however, is coinciden- 
tal. While the separation of the fluid and solid branches 
in the q' 2 diagram is due to the fact that only ordered 
clusters have vanishing q' 2 , there is a large number of 
possible disordered clusters that have q' 6 s=s q^ c , in par- 
ticular, perturbed icosahedral bond arrangements. These 
are, however, not present in the data in large numbers 
and thus can be neglecte d 38 ' 73 . If they occurred in sig- 
nificant abundance in the systems, an increase of (q' 6 ) 
would be the consequence. Values of q' 6 close to gg Cp do, 
however, occur even in disordered systems 3 ^. Thus de- 
viations from q 2 = arguably are a better criterion for 
disorder than deviations from q| cc . 

Since both fee and hep have q 2 = 0, they cannot be dis- 
cerned using q 2 alone. The dense {(f) > 0.649) jLS pack- 
ings, for example, consist of a significant fraction of hep 
and fee clusters on a disordered background. Increasing 
packing fraction reduces the amount of disordered config- 
urations, and proportionally, their weight in the (q[) av- 
erages. Consequently, the q' 2 curves tend towards q' 2 = 
as the ordered clusters take over a larger amount of the 
system, while the terminus of the q' e curves reflects an av- 
erage of qg cc and <Zg Cp , weighted with the relative fraction 
of fee and hep domains. 

A separation of all the regimes can not be seen in the 
q\ plot (d), since q' A takes for crystalline (fee and hep) 
phases fixed values which are lying on a strong random 
background from the disordered parts of the system. 



CONCLUSION 

This article has clearly demonstrated that the con- 
ventional bond-orientational order parameters qi, de- 



fined via nearest neighbor bonds, Eq. ([1}, are very 
strongly affected by the choice of neighborhood definition 
(cf. Fig. [5]); this sensitivity is observed both in the quali- 
tative trend and in absolute values. It was shown that for 
disordered systems without crystallization, <jg strongly 
correlates to the average number of nearest neighbors. 
This effect overshadows the actual structural changes in- 
duced by the physics of the system (cf. FigJ3]). This de- 
pendence is a major drawback that needs to be taken 
into account when using qi for the analysis of particu- 
late matter, especially when comparing qi values across 
different studies. 

We have proposed a unique, well-defined and robust 
structure metric q' l: Eq. that avoids the ambigui- 
ties that come with bond network neighborhoods. Ro- 
bustness of the structure metric is achieved by quantify- 
ing the geometry of the Voronoi tesselation. The MSM 
share the same mathematical form with the conventional 
bond-orientational order parameters, but the "bonds" are 
weighted with the associated Voronoi facet area. This 
guarantees, in particular, that the new Minkowski struc- 
ture metrics are continuous as a function of the sphere 
coordinates. For hep, fee and simple cubic lattices, this 
definition reproduces the values of the conventional qi 
(cf. Tab. H| . For super-cooled hard-sphere fluids, it yields 
a behavior very similar to the conventional q e with the 
(rarely used) rif = 12 neighborhood definition. 

The morphometric neighborhood has previously been 
characterized using Minkowski tensor o 38 ' 53 ' 61 , which 
measure the distribution of normal vectors of the Voronoi 
cells. The Minkowski structure metrics presented here 
can be interpreted as the rotational invariants of a mul- 
tipole expansion of the same distribution of normal vec- 
tors; indeed, the approaches of higher-rank Minkowski 
tensors and Minkowski structure metrics turn out to be 
mathematically equivalent ways to cure the shortcomings 
of bond-orientational order parameters. There are fur- 
ther possibilities to address this problem by introducing 
weighting factors, see for example Ref.~. Note however 
that these approaches need adjustable parameters. The 
caution for the use of qe as a sole determinant of local 
crystallinity expressed in Ref<2£, however, is independent 
of the issues addressed by this paper, and remains valid 
also for the Minkowski structure metric q' 6 . 

Thus, Minkowski tensors and structure metrics both 
provide a "geometrization" of the bond-orientational or- 
der for spherical particles. This suggests a strategy to 
generalize bond-orientational order parameters towards 
aspherical particles, such as ellipsoids, using generalized 
Voronoi tesselations and the q[. Even applications to 
non-cellular shapes with arbitrary topology are possible, 
albeit with altered interpretatio n 64 ! 74 . 

Finally, our analysis supports the more frequent use 
of the low- weight qi, in particular q' 2 , that have been 
largely overlooked in the literature. q 2 carries the same 
information as the anisotropy index f^ 2 of Refsi 38 ' 53 i 61 
(cf. Fig.O. Both q 2 and ft®' 2 can be used to robustly clas- 
sify collective states in particulate matter according to 
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their structural features. Furthermore, q' 2 is very strongly 
discerns between disordered configurations and such of 
high symmetry, such as hep, fee, bec, simple cubic, and 
icosahedral order. 

Clearly, 30 years after the seminal publication by Stein- 
hardt et al, the need for quantitative local structure anal- 
ysis is more evident than ever. The present paper reaf- 
firms the validity and usefulness of the multipole expan- 
sion method. We have, however, described an amended 
version of the bond-orientational order parameters that 
not only renders this method robust and uniquely de- 
fined, but also gives a firmer interpretation of their geo- 
metric meaning. 
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